A glass transition scenario based on heterogeneities and entropy barriers 
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We propose a scenario for the glass transition based on the cooperative nature of nucleation 
processes and entropic effects. The main point is the relation between the off-equilibrium energy 
dissipation and nucleation processes in off-equilibrium supercooled liquids which leads to a natural 
definition of the complexity. From the absence of coarsening growth we can derive an entropy based 
fluctuation formula which relates the free energy dissipation rate in the glass with the nucleation 
rate of the largest cooperative regions. As by-product we obtain a new phenomenological relation 
between the largest relaxation time in the supercooled liquid phase and an effective temperature. 
This differs from the Adam-Gibbs relation in that predicts no divergence of the primary relaxation 
time at the Kauzmann temperature and the existence of a crossover from fragile to strong behavior. 



I. INTRODUCTION 



Standard approaches to the glass transition have been largely based on hydrodynamic or thermodynamic entropic 
theories like that proposed by Adam, Gibbs and Di Marzio (hereafter referred to as AGM) nearly 50 years ago [Q. 
Surprisingly, still nowadays the ideal AGM theory remains not accepted nor disproved, and an unifying description 
of the glass transition problem is still not available. A striking outcome of the AGM theory is the prediction of a 
second-order phase transition driven by the collapse of the configurational entropy (also called complexity) at the 
Kautzman temperature Tr-. 

Strictly speaking an unambiguous definition of the complexity is only possible in mean-field theories where phase 
space splits into ergodic components of infinite lifetime. Recent approaches to the glass transition problem from 
the perspective of disordered systems have, to a large extent, validated this mean-field scenario Q. Nevertheless, a 
complete understanding of the glass transition must go beyond mean field by including nucleation processes into any 
valuable theory. This poses the question about how the mean-field picture for the glass transition is modified in the 
presence of real-space effects. 

Here we propose a scenario for the glass transition where spatial effects due to heterogeneities play the central role. 



II. PHENOMENOLOGICAL DESCRIPTION OF HETEROGENEITIES IN GLASSES 

One of the most intriguing features in glasses is the existence of heterogeneous structures Experimentally 
these manifest as some set of atoms which have a dynamics manifestly slow when compared to the rest. Although 
experiments or numerical simulations in this field are very recent, heterogeneities are a direct manifestation of the 
cooperative nature of nucleation processes Q]. The heterogeneities observed in the experiments are transient frozen 
structures which nucleate on time scales of the order or larger than the observational time. When a droplet of liquid 
nucleates it changes to a new locally disordered structure. The local structure of the glass is always that of a liquid 
and no coarsening of a given pattern is observed. 

Our starting point is tha asspumption that nucleation processes take place everywhere inside the glass when some 
correlated structures of typical size s are built by an anchorage or cooperative mechanism. By anchorage mechanism 
we mean that s atoms must coherently move to occupy positions which enable that region to release some stress 
energy. Each of these moves constitutes an elementary activated process, for instance, the exchange between two 
neighboring particles. In our scenario the time to anchor a region of size s is, 

t(s) /t*\' fBs\ 

where r* = tq exp(i?/T) is the activated time to anchor one atom, tq being a microscopic time and B the corresponding 
energy barrier. 

Experimentally it is known that time correlations in the glass state are stretched but still decay faster than any 
power law. Because of (0) this means that the distribution Ug (t) , the number of cooperative regions (domains) of size 
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s which have already nucleated at time t, must fall down abruptly beyond a maximum size s*{t) for extremely long 
relaxations to occur with a negligible probability. Consequently, ns{t) must have a well defined cuttof at s* such that 
ns{t) for s > s* is nearly zero It can be shown that the existence of this cutoff is tightly related to the cooperative 
character of the dynamics itself eq. and can be illustrated within a simple domain aggregation model |^ . 

Consequently, relaxation to equilibrium of, say, a supercooled liquid is driven by the growth of the largest domain 
size or cutoff value s*{t) since the release of energy to the thermal bath occurs only when the largest domains of size 
s* nucleate for the first time. Nucleation processes of regions of size smaller than s*{t) always occur but do not yield 
a net thermal current flow to the thermal bath. Metastable equilibrium is reached when the cutoff size s* saturates 
and the net energy flow between the glass and the bath vanishes. 

Without knowing any further details about the hydrodynamic equations which describe the relaxation we can set up 
a phenomcnological equation for the energy dissipation inside the glass. This may be described by any thermodynamic 
potential variable, for instance the free energy F. The release of energy per unit time inside the glass is given by 
typical free energy variation of those domains which relax (i.e those with size s*(t)) multiplied by its multiplicity 
Ns*{t) = Vus-it) = (y/s*)n(l) 0, divided by their typical nucleation time t* = exp{Bf3s*), 

1 dF AF*N,.{t) AF* 

V~dt ^ Vt* ^ s*t* ■ ^ ' 

The crucial point which distinguishes nucleations processes in glasses from other physical systems is the fact that 
nucleation occurs between liquid drops. This means that contrarily to standard nucleation theory, the typical release 
of energy A_F* by nucleating a liquid droplet of size s* does not scale with its surface (like in coarsening systems) nor 
its volume (like in standard liquid-solid nucleations) but is finite and independent of s* ^ yielding, 

1^^-^ (3) 
V dt s*t* ■ ^ ' 

where Ap is the average free energy fluctuation. The size s* grows with time as s*{t) ~ Tlog(t) and stops when 
ns{t = teq) reaches the stationary distribution. In principle the value of the equilibration time teq depends on 
the particular model under consideration. From a microscopic theory this information is only contained in the 
hydrodynamic equations. We will see below how we can circumvent the hydrodynamic description by relating the 
relaxation time of the supercooled liquid to the complexity density. 



III. THE FLUCTUATION FORMULA. 



One of the main problems of the standard AGM theory is how to define the complexity Sc- In the original AGM 
theory, Sc was defined as the number of conformational states of the liquid compatible with a given energy. The 
problem in this definition is the following: if the complexity is a static quantity then one is naturally led to the 
conclusion that it vanishes in the supercooled liquid line if measured over time scales much larger that the largest 
decorrelation time. In turn, according to the Adam-Gibbs relation, this implies a divergence of the relaxation time. 
This contradiction can be avoided defining a dynamical complexity Sc{t*, F) by counting all conformations (or basins) 
of the liquid with free energy F which nucleate in a time of the order of the primary relaxation time. Recently, this 
definition has been applied to some mean-field spin glass models where nucleation processes are neglected ||,^. Let 
us assume that the liquid is in off-equilibrium after a quench to temperature Tf. Here we define Sc(t*,F) as the 
logarithm of the total number of conformations with free energy F which can release energy after nucleating for 
their first time at time t* . In other words, Sc{t*,F) counts the number of still not visited conformations at time t*. 
In this scenario regions which have nucleated once have already released their energy while regions which have not 
nucleated still contain some strain energy. Nucleations which do not lead to a release of strain energy do not yield 
new conformational states. 

Each of these conformations at time t* may contain different possible configurations C which do not contribute to 
the complexity but contribute to the free energy of the conformation B, FB{t*,Tf) = — T/ log(^ggg exp(— i?(C)/T/) 
allowing also to define time dependent expectation values of any observable Ai3{t*,Tf) as follows: Ai3{t*,Tf) = 

Ec,«^(C)exp(-m.^(^). 

The liquid character of the glass phase implies that basins with identical free energy FB{t*,Tf) have the same 
physical properties. Note that the free energy and the expectation value of any observable A evaluated at a given con- 
formation B depend on both Tf and the quenching time i*through the configurations C contained in B In particular, 
note that for large times basins with very high free energy only contain configurations which have nucleated several 
times, hence they do not contribute to the complexity. Contrarily, conformations with very low free energy contain 
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configurations which have still not nucleated so they contribute to the complexity. This leads to a dynamical complex- 
ity having a time-dependent cutoff value F* such that the number of conformations or basins with Fsit* ,Tf) > F* 
vanishes. The complexity is a monotonous increasing function of the free energy since the number of possible confor- 
mations which decrease the energy is larger when nucleating regions are smaller. In the asymptotic large t* limit it 
takes the form, 

S,{t*,F< F*) = S,{F,Tf) ; S,{t*,F> F*) ^ -oo (4) 

Since at time t* after the quenching nucleations inside the glass occur between cooperative regions with a local 
liquid structure and no coarsening process drives the nucleation, nucleations must be entropically driven. In other 
words, the probability to jump from any conformation with free energy F onto a conformation with free energy F' is 
always proportional to the number of configurations with final free energy F' , 

WF^F"x^l^^eMSc{F',Tf)-S,{F,Tf)) . (5) 
Substituting in (||) and denoting SF = F' — F we obtain for the transition probability, 

WsF oc exp[/3eff(F*)^i^j e{-6F) (6) 

where we have defined the effective temperature l3csiF*) — as the density of complexity with free 

energy F* . This is the fluctuation formula which establishes the probability of free energy fluctuations. Fluctuations 
to conformations or basins which increase the free energy are forbidden, simply because they have already nucleated in 
the past. While fluctuations to very low free energy conformations are entropically suppressed due to the monotonous 
increasing property of Sc{F, Tf). 

Using the above expression for the free energy fluctuations it follows that in the off-equilibrium state the average 
rate variation of the free energy at time t* is then given by, 

IdF J^^xW^dx _ 1 

V dt ^ t*J^^W,dx /3eff(F*)i*- ^' 

Consistency between this equation with the hydrodynamic equation (|^) requires that s* = Ap(3cs{F*)i i-e. the size 
of the largest nucleating regions s* is directly proportional to the complexity density (also called effective temperature) 
f3es{F*) evaluated at the time-dependent free energy F* . 

From the scenario drawn above off-equilibrium relaxation goes as follows: as time goes on, F* decreases up to time 
tcq when it converges to the equilibrium value Foq — F{toq,Tf). At this time, s* has saturated to a value s*{teq) 
which determines the phcnomcnological relation, 

r„ = To exp(Bs*(teq)/T/) = tq exp{BAp(3,s{F,q)/Tf). (8) 

In the AGM scenario [Q the number of different conformations at time t* corresponds to all possible combinations 
obtained from the two possibilities (nucleated or not nucleated) for the different independent largest nucleating 
regions ns-*{t*). The number of conformations is given by fij. = 2^ yielding for the complexity, Sc{t*) = log(rif) = 
V{\og{2)/s*) directly relating the size of the cooperative region to the complexity. Using eq.(^ and taking t* = t^q 
this yields the famous Adams-Gibbs relation Tq, — exp{B\og{2)/TfSc). By analyzing in details the AGM scenario we 
see that it contains the strong assumption that in the off-equilibrium regime the glass explores conformations differing 
by the nucleation of the largest regions as if they were in equilibrium at the quenching temperature Tf, assumption 
which obviously cannot hold if the system is in off-equilibrium and free energy fluctuations are biased towards lower 
free energy conformations. Note that according to AGM, the size of the cooperative region diverges at the Kauzmann 
temperature Tk while in the scenario presented here it saturates to a flnite value. 



IV. FLUCTUATION-DISSIPATION RELATION 

One more consequence of the fluctuation formula (^ concerns the fluctuation-dissipation ratio (FDR) and its one- 
step character [^O[0 . After quenching to Tf a possible way to quantify violation of FDR is to measure the average 
value of any observable A after a perturbation field Ha conjugated to the observable A is applied to the system. Due 
to the liquid order of the cooperative regions, if the perturbation does not bias one class of conformation over others. 
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the entropically driven assumption implies that states with free energy F are sampled with a probability according 
to their number. In the presence of a field Ha the complexity must be a function of three variables, Sc{F, Tf, A). We 
can simply obtain the average change in the expectation value < A{t) > after switching on the perturbation field Ha 
at time ■ In the linear response regime the Onsager regression principle implies for the transition probabilities 

mF,SA=A-Ao - mp eMiSciF, Tf,A)- S,{F, T^, Aq)] 

.W„e.„(i«^) , (9) 

V ' Ao 

Using the relation (^^^x^)^o = —Pes{F)hA and expanding for Ka small we finally get the famous violation FDT 
expression, 

= Mn{{A{t){A{t) - A(i^))),,=o. (10) 

The description of the violation of FDT in terms of a single time scale t* is therefore consequence of the asymmetric 
shape (i.e. Wsf>o = 0) of (i). 



V. MAIN IMPLICATIONS OF THE PRESENT SCENARIO. 



It can be proven that pcs{F) < 1/Tk where Sc{Tk) = at the Kauzmann temperature Tk- This implies an 
asymptotic crossover for all fragile liquids to strong behavior. A strong glass is a fragile one which has exhausted all its 
complexity and the effective barrier has saturated to its maximum value /3cff = Instead fragile glasses have high 

excess complexity and still big variation of the complexity density f3cS along the supercooled line. The dependence 
of the activation barrier in (^) through the non-universal quantity /3cff explains why it is so difhcult to find a unique 
empirical law that properly describes the viscosity anomaly of all glasses. The prediction that s* saturates to a finite 
value s* = Ap/Tk at Tk is not easy to check experimentally due to the difficulty to cover one order of magnitude 
in the effective barrier. Despite that experimental results suggest a crossover from fragile to strong crossover scenario 
the present scenario can be also checked doing numerical aging experiments. According to relation (|^) the 
t^, dependent effective temperature Pcsitw) (measured through FDT- violations or the formula ^) should be given 
by Pcs(tw) ~* (Tf /BAp)log{tiij/To) if Taitw) — tw and both tq and BAp are nearly and Tf independent. This 
relation should hold for all Tf and predicting a value for the activation barrier and the cooperative size ]T6| . The 
experimental confirmation of a saturation of the heterogeneity sizes for not too fragile glasses and the crossover from 
fragile to strong behavior is probably not out of reach and would be a check of the validity of the present theory. 
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